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Context. Recent measurement of mass of PSR J1614-2230 rules out most of existing models of equation of state (EOS) 
of dense matter with high-density softening due to hyperonization, based on the recent hyperon-nucleon and hyperon- 
hyperon interactions, leading to a "hyperon puzzle". 

Aims. We study a specific solution of "hyperon puzzle", consisting in replacing a too soft hyperon core by a sufficiently 
stiff quark core. In terms of the quark structure of the matter, one replaces a strangeness carrying baryon phase of 
confined quark triplets, some of them involving s quarks, by a quark plasma of deconfined u, d, and s quarks. 
Methods. We construct an analytic approximation fitting very well modern EOSs of 2SC and CFL color superconducting 
phases of quark matter. Then, we use it to generate a continuum of EOSs of quark matter. This allows us for simulating 
continua of sequences of first-order phase transitions, at prescribed pressures, from hadronic matter to the 2SC, and 
then to the CFL state of color superconducting quark matter. 

Results. We obtain constraints in the parameter space of the EOS of superconducting quark cores, EOS.Q, resulting 
from M max > 2 M©. These constraints depend on the assumed EOS of baryon phase, EOS.B. We also derive constraints 
that would result from significantly higher measured masses. For 2.4 Mq required stiffness of the CFL quark core should 
have been close to the causality limit, the density jump at the phase transition being very small. 

Conclusions. Condition M max > 2 Mq puts strong constraints on the EOSs of the 2SC and CFL phases of quark matter. 
Density jumps at the phase transitions have to be sufficiently small and sound speeds in quark matter - sufficiently 
large. A strict condition of thermodynamic stability of quark phase results in the maximum mass of hybrid stars 
similar to that of purely baryon stars. Therefore, to get M max > 2 Mq for stable hybrid stars, both sufficiently strong 
additional hyperon repulsion at high density baryon matter and a sufficiently stiff EOS of quark matter would be 
needed. However, it is likely that the high density instability of quark matter (reconfinement of quark matter) indicates 
actually the inadequacy of the point-particle model of baryons in dense matter at p > 5 ~ 8po. 
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1. Introduction 

The mass of PSR J1614-2230, 1.97 ± 0.04 M 
(jDemorest et all I2010D . puts a constraint on the equation 
of state (EOS) of dense matter in neutron star (NS) cores. 
Namely, maximum allowable mass calculated using an 
acceptable EOS, M max (EOS), should be greater than 
2.0 Mq. This proves crucial importance of strong interac- 
tions in NS cores: their repulsive effect triples the value 
of M max compared to that obtained for non-interacting 
Fermi gas of neutrons, 0.7 M Q . 

Observational constraint M max (EOS) > 2.0 M Q is easy 
to satisfy if neutron star cores contain nucleons only, and 

realistic nuclear forces are used, and we get AfSc(EOS) > 
2.0 M w for m any realistic nucleon interaction models 
(|Lattimerl . l201ll ). However, nuclear interaction models, con- 
sistent with experimental data on hypernuclei, predict the 
presence of hyperons at the densities exceeding 2 — 3/?o, 
where normal nuclear density po = 2.7 x 10 14 g cm~ 3 
(baryon number density uq = 0.16 fm -3 ). Hyperonization 
of the matter implies a softening of the EOS, due to replac- 
ing of most energetic neutrons by massive, slowly moving 
hyperons. For realistic models of baryon interactions one 

< 1.5 Mq (see, e.g. iBurgio et afll201ll 



gets then M^c 



IVidana et al.ll201lT ISchulze fc Riikenll201ll and references 
therein). Such a low Af max is only marginally consistent 
with 1.44 Mq of the Hulse- Taylor pulsar, but was con- 
tradicted already by 1 .67 ± 0.04 M Q of PSR J1903-0327 
([Champion et alJ2 008 ; more p recise value has been recently 
obtained by Freire et al.ll201ll ). 
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Two solutions of the problem of a too low Mm ax have 
been proposed after the discovery of a 2 Mq pulsar. 

Strong hyperon repulsion at high density. It has been 
suggested that adding a new component to the hyperon- 
hyperon interaction, important for p > 5po, can stiffen the 
high-density EOS.B sufficiently to yield M max > 2.0 M Q . 
Repulsive interaction between baryons is supplied by the 
exchange of vector mesons (spinal). Hyperon repulsion due 
to exchange of vector (f> mesons allows for M max > 2.0 Mq, 
without spoili ng the agreemen t with nuclear and hyper- 
nuclea r data (|Bednarek et all I2012L IWeissenborn et all 
l2012afllLastowiecki et al. ll2012t ). lDexheimer fc Schramm! 
(|2008l ) give an earlier general discussion of vector-meson 
contribution to EOS. An additional increase of M max 
(above 2.1 — 2.2 Mq) can be obtained via some break- 
ing of the SU(6) symmetry, usually applied to generate 
vector-mes on - hyperon coup l ing con stants from the nu- 
cleon one (| Weissenborn et all . l2012bl ). Hyperon repulsion 
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at high density is limited by the con dition of thermody- 
namic stability (jBednarek et all |2012() . 

Stiff quark cores in NS. From the point of view of 
quantum chromodynamics (QCD), appearance of hyper- 
ons in dense matter is associated with presence of the 
s-quarks, in addition to the u and d ones confined into 
nucleons. Some authors suggested that the the hyperon 
core in NS could ac tually be replaced by a core of the 
u-d-s quark matter (|Baldo et al.l 120061 and references in 
ISchulze fc Riikenl 12011ft . Let us denote the EOS of the u- 
d-s quark matter by EOS.Q. To yield M max > 2.0 M , 
quark matter should have two important (necessary) fea- 
tures: (1) a significant overall quark repulsion resulting 
in a stiff EOS.Q; (2) a strong attraction in a particu- 
lar channel resulting in a strong color superconductiv- 
ity, needed to make the deconfined Q-phase energetically 
preferred over the confined B(baryon) one. After the an- 
nouncement of the discovery of a 2 M© pulsar, sev- 
eral models of quark cores of NS (hybrid stars), hav- 
ing the properties necessary to yield M max > 2.0 M©, 
have b een proposed ( Ozel et all I2010L IWeissenborn et al.l 
20111 iKlaehn et all |2011L iBonanno fc Sedrakianl 1201 



Lastowiecki et al. 112012ft . The EOS of the hybrid baryon- 



quark (BQ) stars (EOS.BQ) was constructed using a two- 
phase model of B — >Q transition, with different underlining 
theories of the B and Q phases. 

There exist many mode ls of color superc onducting 
quark matter states (see, e.g.. lAlford et al.ll2008ft . Two ba- 
sic states are: two-flavor color superconducting (2SC) state 
and color flavor locked (CFL) superconducting state. In the 
2SC state only light u and d quarks are paired. The 2SC 
state is predicted to be the ground state of quark matter 
at rib 4n . On the other hand, the CFL superconductor, 
in which all three flavors are paired, is predicted to prevail 
at high density, <; Anp. Oth e r sup erconducting states 
are also predicted (jAlford et al.l . 12008ft . but they will not 
be considered here. 

In the present paper we derive constraints on the 
EOS.BQ using an analytical description of EOS.Q. This 
allows us to consider a continuum of the EOS.Q models, 
or, for a given EOS.B, a continuum of the EOS.BQ mod- 
els. As the other authors, we use a two-phase description of 
1st order phase transitions (no mixed-phase state;, neglect- 
ing possibility of a mixed-phase s tate does not influ ence 
much the value of M max , see, e.g., I Afford et al.ll2005ft . We 
also discuss thermodynamic stability of the Q phase in stiff 
quark cores and its i mpact on M mal (this prob lem has been 
already mentioned in lLastowiecki et al. II2012T) . 

EOSs of baryon matter, used in our work, are de- 
scribed in Sect[2J In Scct|3] we construct an analytic ap- 
proximation of the EOS of quark matter, and we show that 
it gives a precise approximation of several existing models of 
color-superconducting quark cores in neutron stars. Our an- 
alytic approximation is then used to construct a continuum 
of EOS.Q, suitable for constructing an EOS.BQ with phase 
transitions at prescribed pressures, for a given EOS.B. In 
Sect 13. 21 we construct continua of the EOS.BQ models, for 
two models of EOS.B: a soft one, significantly violating 
a 2.0 M© bound, and a stiff one, satisfying this bound. 

Constraints on EOS.Q, resulting from mS' > 2 M©, are 
derived in SectJU We then study, in ScctfSl the thermody- 
namic stability of a stiff high-density quark matter. Finally, 



in Sect|5]we summarize and discuss our results, and point 
out the weak points of our models. 

Preliminary results of our work were presented at the 
ERPM Pulsar Conference, Zielona Gora, Poland, 24th-27th 
April, 2012. 

2. EOSs of baryon matter 

There are two types of existing EOS.B, leading to M max < 
2 M© and M max > 2 M©, respectively. They will be here- 
after referred to as soft baryon EOS and stiff baryon EOS. 
We will select one EOS belonging to each of these two 
groups, and use them to illustrate the features characteris- 
tic of the two types of EOS.B. 



2.1. Soft baryon EOS 

Most of existing EOS.B yield A/max significantly lower than 
2 M©. For them replacing a soft hyperon core by a stiff 
quark one seems to be an only way of getting M max > 
2 MfrEI As an example, w e consider a very recent EOS.B 
of ISchulze fc Riikenl (|2011ft . This EOS was obtained using 
Brucckncr-Hartree-Fock many-body approach, and a real- 
istic up-to-date baryon interaction. In the nucleon sector , 
Argq nne Vis nucleon-nucleon potential (|Wiringa et al. I . 
[1995ft was used, s upplemen t ed w ith a phenomenological 
three-body force (|Li et al. I . l2008|) . The hyperon- nucleon 
and hyperon-hyperon potentials were Nijmegen ESC08 

dRiiken et al.ll2010allbh. Thi s EOS.B gives a low Af^ax = 
1.35 M^ (ISchulze fc Riikenl 12011ft . It will be hereafter re- 
ferred to as SR one. 



2.2. Stiff baryon EOS 

The number of EOS.B satisfying Af max > 2 M© started 
to increase steadily after the disc overy of a 2 M ffi pulsar . 
We selected the BM165 EOS.B of lBednarek et all (|2012ft . 
This EOS was obtained using non-linear relativistic mean 
field model involving the baryon octet coupled to meson 
fields. The effective lagrangian includes, in addition to 
scalar and vector-meson terms, also terms involving hidden- 
strangeness scalar and vector-meson coupled to hyperons 
only. For this EOS.B we get M^X = 2.04 M©. It will be 
hereafter referred to as BM165. 



3. Quark-matter cores and their EOS 

We consider electrically neutral u-d-s quark matter in beta 
equilibrium and at T = 0. The baryon number density 



3 ( n u 



n d 



3 ), the energy density is denoted £, 
matter density p — S/c 2 , and baryon chemical potential 
fih = dS/dn^. An important relation is rib = dP/dph- The 
phase transitions under consideration are assumed to be 
1st order ones. Therefore, they occur at a specific (sharp) 
value of pressure, and are accompanied by a density jump, 
from pi to p2- This is a good approximation for B — >Q 
phase transition in cold dense matter, because the smooth- 
ing effec t of the mixed B- Q state is small and can be ne- 
glected (jEndo et al.ll2006l and references therein). An im- 



1 It seems that such a possibility was first considered in 
(|Baldo et al. I. [20061 s ): we are grateful to David Blaschke for call- 
ing our attention to this paper. 
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portant parameter characterizing the density jump at the 
interface between the two phases is A = pil p\. Actually, as 
we have already mentioned in SectflJ inclusion of a mixed- 
phase state does no t change, to a very g ood approximation, 
the value of M max (jAlford et all . 120051) . 
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Fig. 1. (Online color) Calculated points of EOS o f color su- 
perco nducting quark matter in the P — pb plane (|Agrawall 
12010(1 and their approximation by our analytical formula. 
See the text for additional details. 



200 




1000 



1100 



1200 
M„ [MeV] 



1300 



1400 



Fig. 2. (Online color) Calculated points of EOS of color 
superconducting qua rk matter in the P — pb plane 
(jBlaschke et all [2010(1 and their approximation by our an- 
alytical formula. See the text for additional details. 



3.1. Analytical approximation - 25C and CFL phases 

Our method is based on observation that starting from a 
simple linear formula one is able to get a rather precise 
analytic representation of modern EOS.Q(iS) in a phase 
>S =2SC,CFL, of color superconducting quark matter under 



conditions prevailing in neutron star cores. In what follows, 
we will usually omit, for simplicity, the phase label S. 

Linear EOS {P being a linear function of p) is char- 
acteristic of a simplest bag model of quark matter that 
assumes massless quarks, but it also holds with very high 
accuracy for m ore realistic bag model with massive s-quark 
(IZduniM . l2000h . 

Linear EOS is determined by three parameters: a. £ t 
and n* , where a is square of sound velocity in the units of 
c, and £ * and n* are energy and baryon number density at 
zero pressure, respectively. We get then: 



P{£) 


= a(£- 


-£.) 




= p* ■ 




n(P) 


= n* • 


i + i±^f 



0/(1+0) 

1/(1+0) 



l/o 



(1) 



where p* = £ tf jn* is baryon chemical potential at zero 
pressure. The stiffness of the matter is described by the 
parameter a = dP/d£ = (^sound/c) 2 ■ Special cases of linear 
E OSs with a = 1 and a = 1/3 were recently considered 
by IChamel et al. I (|2012f l in their study of exotic cores in 
neutron stars. 

Numerical results for a quark matte r EOS are usu- 
ally given as points in the P — pb plane (jAgrawall . I2010L 
iBlaschke et all 120 ldt) . These varables are very convenient 
to study microscopic stability of matter and for determina- 
tion of phase transition, which correspond to the crossing 
point of P(p) relations for different phases. The density 
jump at phase transition is described then by the change of 
the slope of P(p) function through the relations: n = dPj dp 
and £ = np — P (sec, e.g., Fig 0) 

We assume that the linear EOS, Eq.([T]), accurately de- 
scribes quark matter cores in neutron stars, corresponding 
to baryon density range 2n ^2 rib 10«o Let us stress, 
that the use of linear dependence, Eq.([T]), is restricted to 
2n - 10?7. (or 300 MeV fm~ 3 < £ < 1500 MeV fm" 3 , or 
30 MeV fin -3 < P < 300 MeV fin" 3 ), and by no means is 
claimed to be valid outside the neutron-star core regime. 

Let us introduce dimensionless quantities: rib = nb/n* 
and ~p h = pb/p*- We can then use Eq.([T|) to get P as a 
function of ~p h , 

P=^(^-l), (2) 
v 

where is = (I + a) / a. Then an analytical approximation for 
rib reads 

_ dP(p b ) 



dp h 



Mb 



(3) 



Determination of the value of n* deserves an addi- 
tional comment. It can be taken from original numerical 
calculations if available. In not directly available, it can 
be calculated from the orig inal plot of P< calc '(/j b ) using 
n* = (dP^/dpb)^. 

We now pass to specific cases of <S=2SC, CFL. The 
least-squares fit method results in curves presented in Figs. 
[TJ [2J This fit works very well and could be also checked 
by comparing v alues of n = d P/dp with exact results (if 
available, like in lAg rawal 2 0101) . 

2SC In view of its intermediate-density range, 2no ^ rib ^ 
4no, the 2SC state is less important for M max than the 
high-density CFL state realized for rib > 4n . However, 
as we will show, the softening due to the density jump 
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at the B-Q(2SC) transition has a significant indirect effect 
on the 2 M constraint imposed on the EOS of the CFL 
phase. We consid e red tw o numerica l EOS.Q(2SC)s , calcu - 
lated bv lAgrawall (|2010l ) and two of IBlaschke et all (|2010D . 
All these EOSs were calculated using the Nambu - Jona- 
Lasinio (NJL) model of quark matter and color supercon- 
ductivity. The NJL model is a non-perturbative low-energy 
approximation to QCD. As seen in Figures [TJ [2J our ana- 
lytic formulae fit numerical results very precisely. It should 
be stressed, that these analytical formulae reproduce also 
very well numerically calculated points in t he P — rib plane , 
whenever these points are available, e.g., in lAgrawal (|2010|) . 
Our approximation in this case gives the value of parameter 
v = 4.1 -7- 4.6 which corresponds to a = 0.25 -j- 0.33. 

CFL The baryon density interval Auq ^ rib Ss lOno is cru- 
cial for the value of M max . Therefore, it is the EOS in the 

CFL state which is decisive for the value of Mmejt ■ We con- 
sidered five numerical EOS of CFL superco nducting quark 
matter, t hree models from (|Agrawa 1 l2010Tl and two mod- 
els from (IBlaschke et all l2010f ). All of them were based on 
the NJL model. As we see in Figures [U [5] our analytical 
formulae are very precise. Similarly as for the 2SC phase, 
these formulae reproduce also very well numerically calcu- 
lated points in the P — rib plane. The CFL phase is stiffcr 
than the 2SC one: the values of a range within 0.3 and 0.4 
(y = 3.5 -T- 4.3). 

3.2. A family of analytical models of EOS.BQ 

We generalize now discrete sets {£* ,0^,71**} into a con- 
tinuum of three-parameter models {£3*, a s , n s *}, within a 
region of parameter space determined by appropriate con- 
straints on these parameters. We assume that 0.2 < a < 
0.8, while 2n < n < 10n and 10 MeV fm" 3 < P < 
300 MeV fm- 3 

Having constructed a continuum of the EOS.Q(<S) 
models we are able to simulate, for a given EOS.B, a se- 
quence of phase transition B — >Q(2SC) — ^Q(CFL). 

Transition from B to Q(2SC) Assume that 
B — >Q(2SC) takes place at P = P\. Three parame- 
ters of an EOS.Q(2SC) taken from our family are then 
interrelated by two conditions at the B — !>Q(2SC) phase 
transition point : continuity of the baryon chemical 
potential, and continuity of the pressure, 

M b B) (^i) = vl 2SC) (Pi) , P {B) (£i) = P (2SC) (£ 2 ) • (4) 

Upper indices (B) and (2SC) refer to the baryon phase, and 
the 2SC quark phase, respectively. Now, let us fix the EOS 
of baryon matter, EOS.B. Models of the phase transition 
to the 2SC quark matter are then labeled by Pi = P^(£i) 
and the relative density jump A 2SC = £2/81, with thermo- 
dynamical parameters (ni, ri2, £1, £2) satisfying Equations 
(fj|. Here, the index 1 refers to the B phase, and 2 to the 
2SC phase. 

Transition from Q(2SC) to Q(CFL) Let us now choose 
the pressure at which the 2SC — >CFL transition occurs, 
P = P2. Using conditions of continuity of the pressure 
and of the baryon chemical potential at P = P2, we ob- 
tain a one-parameter family {EOS.Q(CFL)} attached to 
a specific EOS.Q(2SC) from the previously constructed 
{EOS.BQ(2SC)} family. A continuous parameter within 



{EOS.Q(CFL)} can be a CFL or A CFL . This completes the 
second step in the procedure of constructing a general fam- 
ily {EOS.BQ} with B — >Q(2SC) — >Q(CFL) phase transi- 
tions at prescribed pressures Pi and P2, respectively. 

In order to have a " reference one-phase quark core" we will 
also we consider EOS.BQ(CFL), with a phase transition 
from B directly to the CFL phase of quark matter, at pres- 
sure -P CFL and density p CFL . 

Direct transition from B to Q(CFL) at P CFL . For 
Pi < P CFL < P2 this EOS can be either softer or stiffer than 
in the case of B — >2SC — >CFL, depending on the softness 
of the hyperonic EOS. This is illustrated in Fig.[3]where the 
"middle" model of B — >2SC transition (with A 2SC = 1.15) 
gives a mean stiffness similar to that of the hyperon (B) 
phase (i.e. pressures and densities at the bottom of CFL 
core are almost the same in both cases). 




300 400 500 600 



€ [MeV fm" 3 ] 

Fig. 3. (Online color) Examples of EOS.BQ from a contin- 
uum {EOS.BQ} emerging from a soft SR EOS. Hadronic 
EOSs are plotted as dashed lines: N - nucleon EOS; B - 
nucleons and hyperons. The phase transition B — >-Q(2SC) 
takes place at Pi = 31 MeV fm- 3 , and Q(2SC) — ^Q(CFL) 
at P 2 = 45 MeV fm" 3 . The 2SC phase has a 2SC = 0.302. 
Three examples of B — >Q(2SC) — >Q(CFL) are shown 
corresponding to following choices of {A 2SC , A CFL , a CFL }: 
{1.05, 1.2, 0.57}, {1.15, 1.1, 0.58}, {1.3, 1.1, 0.7}. In all these 
cases maximum mass is equal to 2M Q , i.e., parameters 
{A 2SC , A CFL , a CFL } lie on the bounding curves in Fig. [71 

Several examples of EOS.BQ constructed following the 
procedure described above are shown in Figures [31 S) 

For two considered EOS.B, i.e. SR and BM165, we 
are using two different choices of transition pressures Pi 
and P2. Here, Pi corresponds to the pressure at which in 
original models hyperons start to appear. For SR EOS Pi ~ 
30 MeV fm" 3 , while for BM165 EOS Pj ~ 50McV fm" 3 . 
Thickness of the 2SC phase layer corresponds to AP = P2 — 
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Fig. 4. (Online color) Examples of EOS.BQ from a con- 
tinuum {EOS.BQ} emerging from a stiff MB165 EOS. 
Hadronic EOSs are plotted as dashed lines: N - nucleon 
EOS (black); B - nuclcons and hyperons (red). The phase 
transition B — >Q(2SC) takes place at Pi = 47 MeV fm" 3 , 
a 2SC = 0.302 and A 2SC = 1.3. EOS.Q(2SC) is plotted 
as a blue line. Two examples of Q(2SC) — >Q(CFL) at 
Pi = 69 MeV fm" 3 are shown, corresponding to a stiffer 
CFL phase (a CFL = 0.6, A 2SC = 1.3) and a softer one 
(a CFL = 0.35, A 2SC = 1.5). The EOS.Q(CFL)s are plotted 
as green solid lines. 

Pi « 20 MeV fm -3 . In both cases pressure Pi corresponds 
to rib ~ 2rio. 

4. Constraints on EOS.Q in the a — A plane 

A 1st order phase transition in the neutron star core affects 
the EOS in two ways. First, it softens it because A > 1 
implies some density range with dP/dp = 0. Second, a new 
phase is either stiffer or softer than the less dense phase, 
the stiffness of the quark phase being determined by a. The 
condition M max > 2 Mq imposes therefore a condition in 
the a — A plane. In what follows, we start with a case of 
a purely CFL quark core, and then we consider the quark 
core composed of an outer 2SC layer and an inner CFL 
core. Our results will be illustrated by examples presented 
in Figs[5][7] As in the preceding sections, we considered two 
EOS of baryon matter, a soft SR EOS, and a stiff BM165 
EOS. 

Pure CFL core We consider EOS.BQ(CFL), with quark 
core edge at P = -P CFL . We calculate a locus of points in 
the a CFL — A CFL plane corresponding to M max = 2 M 
This locus is a line Ag^(a CFL ) such that points (a, A) below 
it generate Q-cores satisfying M max > 2 M Q , while those 
above it violate this condition. We can alternatively call 
this locus line aS^fA™, ). Of course the location of lines 

Illlll V tr L / 

°cfl — A crL depends on the value of pressure P CFL (density 
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Fig. 5. (Color online) Constraints in the a — A plane, re- 
sulting from M max > 2 M©, obtained assuming SR EOS 
of baryon matter, and for purely CFL cores (no 2SC layer) 
starting at p CFL = 2.3po? 4.0po- Densities 2.3po an d 3.0po 
correspond to Pi and Pi relevant for the boundaries of 
the 2SC layer in a 2SC+CFL core. Each line is an up- 
per boundary of the region of (a CFL , A CFL ) consistent with 
-Mmax > 2 M Q . These lines are labeled by the density p CFL 
at which a direct transition B — >Q(CFL) takes place. 



/9 CFL ) at which phase transition to quark core (CFL) occurs. 
This dependence is presented in Fig. [5j 

2SC+CFL core, effect of a 2SC layer: general pro- 
cedure We replace an outer layer of the matter in the B 
phase containing hyperons (or a layer of the CFL core with 
pressures P CFL < P < P2), by a layer of the 2SC phase. The 
effect of this additional layer of the 2SC phase depends on 
the value of P CFL . If P CFL = Pi we replace an outer layer of 
the CFL core by the quark matter in the 2SC phase. This 
will result in a softening of the quark core since 2SC phase 
is thought to be softer than the CFL one (due to a much 
smaller superfluid gap), and also because of an additional 
density jump A 2SC . Even if A 2SC is very close to one, we get 
Amax CFL (a) < Ag^(a). However, the effect of 2SC phase 
on the value of Amax CFL is then rather small. With an in- 
creasing value of A 2SC , the effect of softening gets stronger, 
and an allowed region of (a CFL ,A CFL ) gets smaller. In the 
limit of Pcfl = Pi we replace baryon matter by the 2SC 
phase and the net effect depends on the relative stiffness 
of these two phases (see Fig|3]and discussion at the end of 
Sect. 153)1 . 

Stiff EOS.B, 2SC+CFL core The loci <^(A CFL ) are 
shown in Figj6] Although the mass fraction contained in 
the 2SC layer is small, its effect on the size of the allowed 
(a CFL , A CFL ) region is strong. For A CFL > 1.2, required a CFL 
has to be signif i cantly larger than th e values obtained in 
dBlaschke et all . I2010l lAg rawall |2010|) . Simultaneously, at 
°cfl ~ 0-5 the density jump due to the 2SC — >CFL tran- 
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Fig. 6. (Color online) Constraints in the a — A plane, re- 
sulting from M max > 2 M Q , obtained assuming BM165 
EOS of baryon matter, and quark core starting at Pi = 
47McV fm" 3 . Each line is an upper boundary of the region 
of (a CFL , A CFL ) consistent with M max > 2 M© . Solid lines 
are obtained for quark cores composed of a 2SC layer (Pi < 
P < P 2 ) and a CFL core starting at P 2 = 69MeV fm~ 3 . 
These lines are labeled by the density jump at the B- 
2SC interface A ogr , = 1.05,1.15,1.30. We used a 2SC = 0.3 
(|Agrawall l2010f ). Dashed line is obtained for purely CFL 
cores (no 2SC layer) starting at P 2 . 




CFL 



Fig. 7. (Color online) Constraints in the a — A plane, 
resulting from M max > 2 M©, obtained assuming SR 
EOS of baryon matter, and quark core starting at Pi = 
31 McV fm~ 3 . Each line is an upper boundary of the region 
of (a CFL , A CFL ) consistent with M max > 2 M© . Solid lines 
are obtained for quark cores composed of a 2SC layer (Pi < 
P < P 2 ) and a CFL core starting at P 2 = 45 MeV fm~ 3 . 
These lines are labeled by the density jump at the B- 
2SC interface A 2SC = 1.05,1.15,1.30. We used a 2SC = 0.3 
(|Agrawall |2010[) . Dashed line is obtained for purely CFL 
cores (no 2SC layer) starting at P 2 = 45 MeV fm~ 3 . Open 
circles correspond to EOSs presented in Fig. [3j 



sition is con strained to values significantly below the ones 
obtained in (jBlaschke et all I2010L lAgrawalL l2010h . 

Soft EOS.B, 2SC+CFL core The loci <^(A CFL ) arc 
shown in FigJT] Even for a very low density jump A 2SC = 
1.05, we obtain a CFL > 0.4, which is rather stringent. In 
our case, the result obtained for B — ^Q(CFL) at P CFL = 
P 2 is very similar to that of B — >Q(2SC) — s-Q(CFL) with 
A2SC = 1.15 (the dashed line is very close to the solid one) 



5. Stability of quark cores and M max 

Up to now, we did not consider the (thermodynamical) 
stability of a stiff quark core in a hybrid (BQ) star. A 
stiffening of the EOS is necessarily associated with the in- 
crease of the baryon chemical potential (see an example in 
iBednarek et al.ll2012t) . In particular, it may lead to the ther- 
modynamical instability of the stiff (Q) phase with respect 
to the re-conversion into the (B) one. This instability results 
from the violation, above a certain pressure, of the condi- 
tion fjS®\p) < fi^\p). Assuming a complete thermody- 
namical equilibrium, we are dealing with a first-order phase 
transition back to the (B) phase, th at one can call reconfine- 
ment (cf.. lLastowiecki et al. 1120121 ). A corresponding EOS 
will be denoted EOS.BQ and the M(R) branch based on 
this EOS will be labeled BQ. Examples of the N,B,BQ and 



BQ branches in the M—R plane, obtained for a soft SR EOS 
of baryon matter are presented in Fig|5] The reconversion 
Q — > B strongly limits the size of the quark core in hybrid 

stars and results in the value of Mmax ~ Mma X = 1-35 M© 
(Fig® 

For the BM165 EOS.B we get = 2.04 M©. 

Replacing hyperon cores by stiff quark ones can further 
increase the value of M max . An example is shown in FigO 
where we obtain M$£P = 2.07 M©. However, if complete 
thermodynamic equilibrium is imposed, the Q — >B transi- 
tion back to the B phase takes place and one gets maximum 
allowable mass M^S X rj M^S x = 2.04 M©. 

6. Summary, discussion, and conclusions 

The existence of a 2 M© pulsar is a challenge for neutron 
star models with strangeness carrying cores. Strangeness is 
associated with s quark, either confined into hyperons or 
moving in a (deconfined) quark plasma. 

The threshold density for the appearance of hyper- 
ons, predicted by realistic models of dense matter consistent 
with nuclear and hypernuclear data, is ~ 2pg— 3po- Realistic 
baryon interactions lead to M max , for neutron stars with 
hyperon cores starting at such density, that is significantly 
below M = 2.0 M©. This contradiction can be removed 
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Fig. 8. (Color online) Mass M vs. radius R for non-rotating 
neutron star models calculated using different assumptions 
concerning structure of dense matter for p > 2/9o- Dotted 
segments of M(R) curves: configurations unstable with re- 
spect to small radial perturbations. N - nucleon matter , 
Argonne Vis+TBF nucleon interaction ([Schulze fc Riikenl . 
1 2 1 if ) . B - baryon matter, EOS calculated using hyperon- 
nucleon and hyperon -hyperon potentials Nijmegen ESC08 
(|Riiken et all boiOal fbh. Open circles - configurations of 
maximum allowable mass, with central density given in the 
units of po- BQ - hybrid stars with quark cores, when con- 
dition p^(P) < (P) is not imposed. BQ - hybrid stars 

with stable quark cores, when condition p { h Q) (P) < ^ B) (P) 
is satisfied. Analytical EOS of quark matter in the CFL 
state with p 2 = 2.35po, o C fl = 0-5, A CFL = i.2. Enlarged 
rectangle: Vicinity of the high-density 1st order phase tran- 
sition Q — ^B (reconfinement). 



due a hypothetical strong high-density repulsion acting be- 
tween hyperons. As discussed in several papers, this strong 
high-density repulsion could result from the exchange of a 
vector meson </> coupled only to hyperons. 

Howewer, it has also been considered that massive neu- 
tron stars could actually be hybrid stars with stiff quark- 
matter cores that allow for M > 2 M Q . Strong overall repul- 
sion between quarks should be accompanied by a strong at- 
traction (pairing) in a specific two-quark state, correspond- 
ing to a strong color superconductivity with a superfluid 
gap - 100 MeV. 

In the present paper we performed a general study 
of possibility for hybrid neutron stars with quark cores to 
reach M > 2 M Q . We considered a continuum of param- 
eterized EOS of quark matter, including several existing 
models. This allowed us to consider general case of quark 
cores coexisting with baryonic matter at a prescribed pres- 
sure. We determined necessary features of baryon - quark- 
matter phase transition. First, the density at which first- 
order phase transition to quark phase occurs should be sim- 
ilar to the threshold density for hyperons, ~ 2po — 3po- 




Fig. 9. (Color online) Same notations as in FigO 
but for BM165 EOS of nucleon and baryon matter 
(jBednarek et al.l . 12012} ). Analytical EOS of quark matter 



and EOS.Q(CFL) with p 2 = 2.5p , a c 



0.5, A c 



1.2. 



Second the relative density jump at the baryon-quark mat- 
ter phase transition should be below 30%. Third, quark 
matter should be sufficiently stiff, which can be expressed 
as a condition on the sound speed in quark plasma. 

The measured 2.0 M Q is actually a lower bound to a 

true A/ max = M4 t Ix e) . The upper bound, resulting the con- 
dition of speed of sound less than c combined with our 
confidence in the the oretical nucleon EO S for p < 2po, 
is 3.0 M Q (see, e.g., lHaensel et al.l [20071 and references 

therein). Afi'™' lies therefore between 2 M© and 3 M Q . 
Obviously, neutron star masses higher than 2.0 M Q have to 
be contemplated. The question is how much higher? The 
mass of a "black widow" pulsar could be as high as 2.4 M , 
albeit the present uncertainty is too large for this number to 
be us ed as an observational constraint (see, e.g.. lLattimerl 

lint). 

To discuss possibility of reaching masses significantly 
larger than 2.0 M Q we plotted in Fig.[TU]thc bounding lines 

for A/4°a b x S) = 2.2 M Q and = 2.4 M Q . As we see in 

Fig. [TJB to fulfill condition M r ( ° a b x s) = 2.4 M we have to 
assume very stiff quark matter, quite close to the causality 
limit a = 1. 

The situation becomes even more difficult if we require 
a strict stability of quark cores. As a result of high stiffness 
of quark matter, necessary for Mm„' > 2 M ( and a for- 
tiori for higher lower bounds Mmix ' ) , the quark phase turns 
out to be unstable, beyond some pressure, with respect to 
hadronization. Assuming complete thermodynamic equilib- 
rium, we get very similar M max for the stars with hyperon 
and quark cores. Consequently, transition to quark matter 
could not yield M max > 2 M Q if neutron stars with hyper- 
onic cores had Mii (significantly) below 2 M Q . This is 
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0.4 0.6 0.8 1 



Fig. 10. (Color online) Constraints in the a — A plane, re- 
sulting from 3 different values of maximum measured mass 
Mi° a b x s) = 2.0 M ,2.2 M ,2.4 M Q . We assume BM165 
EOS of baryon matter, and quark core starting at P\ = 
47 McV fm -3 . Each line is an upper boundary of the re- 
gion of (a CFL , A CFL ) consistent with M max > A-ff . Solid 
lines are obtained for purely CFL cores (no 2SC layer). 
Dashed lines are obtained for quark cores composed of a 
2SC layer (Pi < P < P 2 ) and a CFL core starting at 
P2 = 69 MeV fm -3 . These lines are labeled by the density 
jump at th e B-2SC interf ace A 9SC = 1.05,1.15. We used 
a 2SC = 0.3 l|AgrawalL[20Tot ). 



true also for Mmax > 2 M Q . Therefore, provided our pic- 
ture of dense matter is valid, we find that a strong hyperon 
repulsion at high density is mandatory in general. 

The high-density thermodynamic instability of the 
quark phase, and its consequences for M max , should be 
taken with a grain of salt. Our models of dense baryonic 
matter assume point particles. This assumption may be 
expected to break down at p ~ 5 8po. Therefore, the "re- 
confinement" of the quark phase is, in our opinion, likely 
to indicate the inadequacy of point-part icle baryonic phase 
model (see also lLastowiecki et al. 112012]) . 

There is another weak point of the commonly used 
models of quark cores in neutron stars, characteristic also 
of the present paper: this is a two-phase approach, with 
each phase, baryon and quark one, treated using basically 
different descriptions. In principle, both phases and transi- 
tion between them should have been treated using a uni- 
fied approach based on the QCD, so that that the influ- 
ence of the dense medium on the baryon structure and 
baryon interactions are taken into account in a consistent 
way. Such approach is beyond the reach of the present day 
theory of dense matter. However, a phenomenological mod- 
eling of baryon structure in dense matter is possible, e.g., 
within a quark-meson cou pling model (for references, see 
IWhittenburv et al. I l2012f) . More complete description of 



neutron-star quark cores, going beyond the two-phase ap- 
proximation, should hopefully be obtained in the future. 

In this paper we were considering non-rotating con- 
figurations. Pulsar PSR J1614-2230 rotates with frequency 
/ = 1/P = 317 Hz and th e effect for max i mum mass is 
of the order of ~ 0.01M o (jBednarek et all . I2012D . much 
smaller than accuracy of mass determination. However it 
should be noted that for neutron star rotating with max- 
imum observed frequency of 716 Hz the effect of rotation 
would be about five times larger. 
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